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L'invention concerne notamment un procede de separation 
autodidacte (ou aveugle), n'exploitant aucune information a priori svlt les 
sources ou sur le front d'onde, a Tordre 4 et aux ordres superieurs a 4, de P 
sources cyclostationnaires (deterministes ou stochastiques, analogiques ou 
numefiques, a modulations lineaires ou non lineaires) et statistiquement 
independantes, a partir d'une reception a N capteurs (N>2). 

Elle s'applique par exemple dans le domaine des 
radiocommunications, des telecommunications spatiales ou de I'ecoute 
passive de ces liaisons, dans des gammes allant par exemple de la VLF a la 
EHF). 

Elle trouve aussi son application dans des domaines tels que 
Tastronomie, le biomedical, le radar, le traitement de la parole, etc. 

La separation autodidacte de sources et plus particulierement 
I'Analyse en Composantes Independantes (ICA) suscitent actuellement 
beaucoup d'interet. Elles trouvent en effet leur place dans de nombreuses 
applications teiles que les telecommunications, le traitement de la parole ou 
encore le domaine biomedical. ' 

Par exemple, en traitement d'antennes, si un certain nombre de 
sources emises sont receptionnees sur un reseau de capteurs, et si pour 
chaque source Tetalement temporel des canaux associes aux differents 
capteurs est negligeable devant le temps symbole, alors un melange 
instantane des sources emises est observe sur les dits capteurs. 

La separation autodidacte de sources a notamment pour but de 
restituer les sources supposees statistiquement independantes et ce 
uniquement a partir des observations regues par les capteurs, 

Selon Tapplication, il est possible de ne retrouver que le melange 
instantane, c'est-a-dire les vecteurs directeurs des sources, ce qui est le cas 
par exemple de la goniometrie oii le dit melange porte a lui seul toute 
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rinformation necessaire a la localisation angulaire des sources : on parle 
alors d 'identification aveugle de melange. 

Pour d'autres applications telles que la transmission, il est 
necessaire de retrouver les sources emises : on emploie alors I'expression de 
separation ou en encore d'extraction aveugle ou autodidacte de sources. 

Certaines techniques connues de I'art anterieur cherchent a 
decorreler (a I'ordre 2) les signaux, comme on peut I'observer en Analyse 
Factorielle avec I'Analyse en Composantes Principales (PCA). 

LICA, pour sa part, vise a reduire les dependances statistiques 
des signaux egaiement aux ordres superieurs. De ce fait HCA permet 
d'identifier, de maniere autodidacte, le melange instantane et par la meme 
occasion d'extraire les sources emises dont au plus une est supposee 
gaussienne. Ceci n'est actuellement possible qu'en respectant certaines 
hypotheses : le melange bruite des sources doit etre lineaire et par ailleurs 
15 surdetermine (le nombre de sources P doit etre inferieur ou egal au nombre 
de capteurs N). 

Si Comon introduit le premier le concept dlCA et propose une 
solution, COM2 dans la reference [1] (les differentes references sont 
regroupees a la fin de la description), maximisant un contraste base sur les 
20 cumulants d'ordre 4, Cardoso et Souloumiac [2], quant a eux, developpent 
une approche matricielle, plus connue sous le nom de JADE, et donnent 
ainsi naissance a I'algorithme de diagonalisation conjointe. 

Quelques annees plus tard, Hyvarinen et al. presentent la 
methode FastICA, dans un premier temps pour des signaux reels [3], puis 
25 dans le cas complexe [4]. Cette demiere introduit un algorithme 
d'optimisation de contraste baptise algorithme du point-fixe. 

De son cote, Comon propose une solution simple, COMl [5], a 
I'optimisation du contraste presente dans [6]. 
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Bien que tres performantes sous les hypotheses enoncees 
precedemment, ces methodes peuvent etre neanmoins tres fortement 



3 



perturbees par la presence de bruit inconnu, Gaussien ou non, correle 
spatialement et inherent a certaines applications comme les 
radiocommunications HF. 

Par ailleurs, comme il a ete dit plus haut, les methodes 
precedentes ne sont congues que pour traiter des melanges surdetermines 
de sources. Or dans la pratique, tel qu'en radiocommunications, la 
reception de plus de sources que de capteurs n'est pas rare, surtout si la 
largeur de bande de reception est importante. On est alors en presence de 
melanges dit sous-determines {P> N}. 

Plusieurs algorithmes ont d'ores et deja ete developpes dans le 
but de traiter ce type de melanges. Certains abordent le difficile probleme de 
Textraction des sources [7-8] alors que le melange ne peut plus etre inverse 
lineairement, tandis que d'autres traitent de Tidentification de la matrice de 
melange [7] [9-12], 

Les methodes proposees dans [9-1 IJ exploitent uniquement les* 
statistiques d'ordre 4 tandis que celle presentee dans [12] repose sur 
rutilisation de la seconde fonction caracteristique des observations, 
autrement dit sur rutilisation des cumulants non nuls de tout ordre. Quant 
a la methode employee dans [7], elle repose sur la maximisation de 
probabilite conditionnelle, en Toccurrence celle des donnees 
conditionnellement a la matrice de melange. 

Si ces methodes sont performantes, elles souffrent d'inconvenients 
en contexte operationneL 

Ainsi la methode [9] est difficile a mettre en oeuvre et ne garantit 
pas ridentification des vecteurs directeurs de sources de meme kurtosis. 
Les methodes [10] et [11] ne permettent pas quant a elles d'identifier les 
vecteurs directeurs de sources circulaires. La methode [10], baptisee S3C2, 
pour sa part confine rutilisateur dans une configuration de 3 sources et 2 
capteurs, interdisant tout autre scenario. La methode [7] autorise 
ridentification de 4 signaux de parole avec seulement 2 capteurs, toutefois. 
les echantillons observes doivent etre temporellement independants et 



chaque source doit avoir une densite de probabilite eparse. Pour finir, la 
methode [12] n'est applicable qu'au cas de sources reelles, ce qui est tres 
restrictif notamment en comraunications numeriques. En outre, 
ralgorithme depend grandement du nombre de sources, et rien ne prouve 
aujourdTiui qu'une mauvaise estimation de ce parametre ne degrade pas les 
performances de la methode. 

La presente invention offre une nouvelle approche reposant 
notamment sur Texploitation de la totalite ou pratiquement de la totalite de 
rinformation propre aux vecteurs directeurs ap des sources, contenue de 
maniere redondante dans la matrice representative des statistiques 
circulaires d'ordre m = 2q du vecteur des enveloppes complexes des signaux 
en sortie des capteurs. 

L'invention conceme un precede didentification aveugle de 
sources au sein d'un systeme comportant P sources et N capteurs. II est 
caracterise en ce qu'il comporte au morns une etape d 'identification de la 
matrice des vecteurs directeurs des sources a partir de I'information 
propre aux vecteurs directeurs ap des sources contenue de maniere 
redondante dans les statistiques circulaires d'ordre m=2q du vecteur des 
observations regues par les N capteurs. 

Les statistiques circulaires d'ordre m = 2q s'expriment par 
exemple en fonction d'une matrice diagonale de rang plein des 
autocumulants des sources et d'une matrice representant la juxtaposition 
des vecteurs directeurs des sources de la maniere suivante : 

Cm.K^ Aq^m^sAq^ (11) 

OU (^m.s = diag{[C/j^^^;;;'^^^ >--m C^;^; ]) est la matrice diagonale de rang plein des 
autocumulants C^^;;;;^,, d^ordre m^2q des sources, de taille (Px P\, et oii Ac^ 
= [ ai®(9-i) ® cif ... a^ici'i) ® ap ], de taille {N^i^Pj et supposee de rang 
plein, represente la juxtaposition des P vecteurs colonnes [ a/ft^-D ® ]. 

Utilisation du precede a un reseau de communication et/ou pour de la 
goniometrie a partir des vecteurs directeurs identifies. 



L'invention presente notamment les avantages suivants : 

• Elle permet ridentification autodidacte de melanges instantanes aussi 
bien surdetermines (le nombre de sources est inferieur ou egaJ au 
nombre de capteurs) que sous-determine (le nombre de sources est 
superieur au nombre de capteurs) ainsi que Textraction aveugle des 
sources dans le cas surdetermine ; 

• A Tordre m==2q, pair q>2, le procede designe BIOME (Blind Identification 
of Over and underdetermined MixturEs of sources) a la capacite de 
traiter jusqu'a P = N «i - sources a partir d'un reseau a N capteurs 
differents , des lors que les autocumulants d'ordre m des sources sont 
de meme signe; 

• Une application du procede BIOME a Tordre 4, nommee ICAR 
(Independent Component Analysis using Redundancies in the 
quadricovariance), permet ridentification autodidacte de melanges 
instantanes surdetermines (P < N) de sources et Textraction aveugle Me 
celles-ci, d'une maniere s'averant robuste vis-a-vis de la presence d'un 
bruit Gaussien inconnu correle spatialement, des lors que les sources 
ont des kurtosis (autocumulants normalises d'ordre 4) de meme signe ; 

• Une application du procede BIOME a Tordre 6, appelee BIRTH (Blind 
Identification of mixtures of sources using Redundancies in the daTa 
Hexacovariance matrix), permet ridentification autodidacte de melanges 
instantanes aussi bien surdetermines {P < N) que sous-determines (P > 
iV) de sources, ainsi que Textraction aveugle des sources dans le cas 
surdetermine, Le procede BIRTH a la capacite de traiter jusqu'a 
sources a partir d'un reseau a N capteurs differents, des lors que les 
autocumulants d'ordre 6 des sources sont de meme signe. 

D'autres caracteristiques et avantages du procede selon 
Tinvention apparaitront mieux a la lecture de la description qui suit d'un 
exemple de realisation non limitatif annexe des figures qui representent : 

• La figure 1 un schema d*exemple de mise en oeuvre du procede a I'ordre 
m, 

• Les figures 2, 3 ef 4 des resultats de simulation de la mise en oeuvre du 
procede a Tordre 4, 



• Les figures 5, 6, 7 et 8 des resultats de simulation pour le procede 
applique a Tordre 6. 

Les exemples qui suivent sont donnes pour ridentification et/ou 
rextraction de sources dans un reseau comportant une antenna reseau 
comprenant N capteurs. Cette antenne est supposee recevoir un melange 
bruite de P sources par exemple a bande etroite (BE) et statistiquement 
independantes. 

Sous ces hypotheses, le vecteur x{t) des enveloppes complexes 
des signaux en sortie des capteurs s'ecrit, a Tinstant t 

p 

x{t) - X ^P^^^ + v<t) = s(t) + v[t) (1) 
p=l 

• ou v(t) est le vecteur bruit, suppose centre, Gaussien et inconnu, 

• Sp{t) et ap correspondent respectivement a Tenveloppe complexe, BE, 
cyclostationnaire et cycloergodique (avec un eventuel residu de porteuse 
le cas echeant), et au vecteur directeur de la source p, 

• s(t) est le vecteur dont les composantes sont les signaux Sp(t) et A est la 
matrice (N x P) dont les colonnes sont les vecteurs ap. 

Par ailleurs, le procede enonce de maniere generale ci-apres a 
Tordre m 2g (q > 2) utilise les hypotheses suivantes Hi-^: 
HI : A tout instant t, les sources a valeurs complexes Sp(t) sont 
cyclostationnaires, cycloergodiques et mutuellement decorrelees a Tordre 

m ; 

H2 : A tout instant ^, les composantes Vrt(t) du bruit sont stationnaires, 
ergodiques, gaussiennes et circulaires; 

H3 : A tout instant f, s[t) et v(i) sont statistiquement independants ; 

H4 : les autocumulants d'ordre m des sources sont non nuls et de meme 
signe. 



Avec les hypotheses precedentes, pour un ordre pair m = 2q 
donne, le probleme de ridentification autodidacte de melanges instantanes 
de sources consiste a trouver, a partir de rexploitation de certaines 
statistiques d'ordre m des observations, la matrice A a une matrice triviale 
pres (une matrice triviale est de la forme A n ou A est une matrice diagonale 
inversible et n une matrice de permutation). 

La separation (extraction) aveugle de sources, consiste 
notamment a determiner le separateur Lineaire et Invariant dans le Temps 
(LIT), W, de dimension [N k F), dont le vecteur de sortie, de dimension (Px. 
1), 

y(t) - WH x(t) (2) 

correspond, a une matrice triviale pres, a la meilleure estimee, s(/) , du 
vecteur s(t), ou le symbole " signifie transpose conjugue. 

Le separateur W est defini a une matrice triviale pres dans la 
mesure oii ni la valeur des puissances de sortie, ni I'ordre dans lequel les 
sorties sont rangees ne changent la qualite de restitution des sources. 

Avant d'exposer les etapes du precede selon Tinvention, quelques 
rappels sur les statistiques des observations sont donnes. 
Statistiques des observations 

Le precede selon Tinvention utilise notamment les statistiques 
circulaires d'ordre pair, 2, 4, 6, m des observations. 

Selon Tart anterieur decrit dans la reference [15], I'expression des 
cumulants dbrdre m en fonction des moments dbrdre inferieur a m peut 
etre sihaplifiee. ^ 

une, quantite scalaire a valeur complexe dependant des q 
indices inferieurs d, e et des q indices superieurs /, g a valeur dans 
{1, 2, iV}. La quantite G/;;:;;f^ verifie alors les trois S3nnetries suivantes : 
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• toute permutation des indices inferieurs de G^""f^ entre eux ne 
modifie pas la valeur de G^'j'l^ : par exemple, pour q = 2, G//^ = 

• toute permutation des indices superieurs de G^' 'f^ entre eux ne 
modifie pas la valeur de G/;;;^^, : par exemple, pour q = 2, G^'{^- 

• permuter tous les indices superieurs avec tous les indices inferieurs 

^iZtx ^ po^r effet de conjuguer la valeur de G^^;' f^^ : G^' ^ ^ 

D 'autre part, la notation suivante est adoptee : la quantite 
i^^^k^tx designe la combinaison lineaire des r produits 

possibles et distincts (modulo les trois symetries decrites precedemment) de 
type ^j;:;;^ G(;:i^, ... g;:;% , ponderes par la valeur 1. Chacun des r produits 

est construit a partir du produit G^;;::;l^ GI::1^.,,G;'-^^ utilisant et en 
combinant les deux regies de permutation suivantes : 

• un indice inferieur d\in des termes du produit G{ ^ G/ * G," 
permute avec un indice inferieur d'un autre terme du meme produit 
GdZj.K ^hZtx'"G";^^ pour donner un autre produit (distinct) : par 
exemple, pour (qi. qa) = (2, 2), G/;^, G^^^donne comme autres 
produits distincts (modulo les trois symetries decrites precedemment) 

• un indice superieur d'un des termes du produit G/ ^ Gl"'"^ G" 
permute avec un indice superieur d'un autre terme du meme produit 
^iztx Gr,Ztx'-G";;^^ pour donner un autre produit (distinct) : par 



9 



exemple, pour (qi, ^2) = (2, 2), G^;f^G/^f^donne comme autres 
produits distincts (modulo les trois symetries decrites precedemment) 

alin d'obtenir la totalite des r produits possibles et distincts (modulo les 
5 trois symetries decrites precedemment). 

L'exemple suivant, ou (qi, 92) = (2, 1) et r = 9, illustre cette notation : 
+ G-J^e,. Gl^ + G|;;., Gl + Gj,->., Gl^ + G/^^ G* , + G}'* GJ,, (3) • 



10 Statistiques d'ordre m 

Dans le cas de sources stationnaires ou cyclostationnaires 
potentiellement non centrees, les statistiques circulaires d'ordre m du 
vecteur x(f), donne par (1), s'ecrivent 



.15 



OU le symbole < J{t) >c = (1/T) f' f[t)dt correspond a I'operation de 

J-r/2 



rr/2 
-r/2 

moyennage temporel en temps continu de J{t) sur un horizon d'observation 
infini. 



^S^"'"' ^'>=^""^^ <'>*'V 

20 (5) 

ou q termes sont conjugues et q termes sont non conjugues. 

Les statistiques d'ordre m decrites par I'expression (5) sont dites 
circulaires car le cumulant Cum{xa((l, Xe{ti,..., x^Q, Xg{t]*, 
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d'ordre m est calcule a partir d'autant de termes conjugues {xg(t)*, Xh{t)*y 
Xi{t)*) que de termes non conjugues (xd{t), Xe(t), xj{t)), 

Les cumulants circulaires d'ordre m peuvent etre exprimes en 
fonction de moments d'ordre inferieur a m de la maniere donnee ci-apres. 

Soient m'^;^^'V^' "''"(0 et C^+^'V^'' les moments et 

cumulants circulaires d'ordre m associes au vecteur dbbservation x{t), 
definis par : 



^l^t^'H-;^X ''"(^) = E[X. (t),^. (t),...,X. Xt),X. (tr,X. (t)*,...,X. (t)*] 

^«<x = fe,(l)U6';_,(l)U...U5'g,,(/r)} designe la « g ..-ieme partition, parmi 
« Gk » possibles, de « k » sous-ensembles Sa \ =^"!^' ■ 'f^ , oil ^ < r ^ s * ... * t 



< g et Q+1 < u ^ 1/ ?i ... ;t w < m. telle que Par/^ ^ =^t^'^+2' q^^^ ^ 

i 2'" 9 

Toperateur union, note (J , il verifie 



De maniere pratique, deux grandes families d'estimateurs 
peuvent etre utilisees pour estimer les statistiques precedentes : dans le cas 
de sources stationnaires, on peut employer un estimateur empirique non 
biaise et consistant pour des sources stationnaires, ergodiques, alors que 
dans le cas de sources cyclostationnaires potentiellement non centrees, il 
faut employer un estimateur dit exhaustif, non biaise et consistant pour des 
sources cyclostationnaires, cycloergodiques, potentiellement non centrees. 
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Get estimateur exhaustif est par exemple determine selon une demarche 

decrite dans les references [13-14], 

Agencement et stockage des statistiques d'ordre m 

Comme il est montre ci-dessus les statistiques C^' '""'y. , pour 

5 1 < d, e, f; g, h, k < iVsont des fonctions a m = 2q entrees {m pair), il 
alors possible de les ranger dans une matrice [N^ x N^) que Ton nommera 

Cm,x. 

De maniere explicite, la quantite C^'^^' j ^ se trouve situee a la 

eme ligne et a la j-eme coloime de la matrice Cm,x en posant i = iV[...iV[iV(d- 
10 .,.]+fcet j= N[ ,.,NlN{g-l) + hr-l]+ ...]+/ 

Agencement et statistiques d^ordre 4 

Dans le cas de sources stationnaires centrees, les statistiques 
circulaires d'ordre 4.du-vectetir x(t), donne par (1), s^ecrivent 

15 Ci;t = Cum{x4t). xe(t). xAt)\ x,m - M/X, - , , M/'^ ^ 12] M^^ (4)' ^ 

De maniere pratique, ces statistiques peuvent etre estimees en employant 
un estimateur empirique . non biaise et consistant pour des sources" 
stationnaires, ergodiques, centrees. 

20 Dans le cas de sources cyclostationnaires potentiellement non 

centrees, les statistiques circulaires d'ordre 4 du vecteur a prendre en 
consideration, s'ecrivent 

CfrL ^ <.Cum{xe^(t), xe[t], xAt)\ x^m >c = < M/^^ W >c - < \^M:,^i^M{f{t\ 
>c - < , , [t) M(^' {t)>c'< [2] M/, [tiMl (t) >c 2 < ^ (t) A/^ ^ (f) M/'^ (£) 
25 >c 2.< M/{£)Mf U)M,,,(t) >c + 2 < [4]^^^^ (^)M/^, (t) >c - 6 < 

^rf,x it) A/,,, (t) M,^ (t) Ml [t] >c {Sy 
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De maniere pratique, ces statistiques d'ordre 4 peuvent etre 
estimees en employant Testimateur dit exhaustif non biaise et consistant 
pour des sources cyclostationnaires, cycloergodiques, potentiellement non 
centrees. Get estimateur exhaustif est decrit dans [13-14]. 

Comme il est montre ci-dessus les statistiques C//, , pour 1 ^ d, 
e,f,g<N sont des fonctions a quatre entrees, il alors possible de les ranger 
dans une matrice (iV2 x JV2) que I'on nommera matrice de quadricovariance 
Qx donnee par exemple dans la reference [13-14]. 

Agencement et Statistiques d'ordre 6 

Dans le cas de sources stationnaires centrees, les statistiques 
circulaires d'ordre 6 du vecteur x(t), donne par (1), s'ecrivent 

c^tf.. = cum{xd(g, xe{t), xj{t^, xg[r, XH{t\*, xm = ^1;';., - mMl^^^M'/ - 

(6) 

En pratique, ces statistiques peuvent etre estimees en employant 
un estimateur empirique non biaise et consistant pour des sources 
stationnaires, ergodiques, centrees. 

Dans le cas de sources cyclostationnaires centrees, les 
statistiques circulaires d'ordre 6 du vecteur x({), donne par (1), s'ecrivent 

20 Qe!;.x = < Cum{jQ(f), x«(f), xAt^, xg{t)*, xh[t)*. Xi(t)*} >c = < >c - < 

l^WLfJt)M'/[t) >c- < [9JM^;j;jt)M;,, (t) >c- < [3]M,^Jt)Mf-'/ (ti >c+ 2 < 

f 9) ^ (t) M;,, {t)M';'(t)>c+2< [6] Ml^ (t) M {t) M'^^ (^) >c (7) 

De maniere pratique, ces statistiques d'ordre 6 peuvent etre 
estimees en employant un estimateur dit exhaustif non biaise et consistant 
25 pour des sources centrees cyclostationnaires, cycloergodiques. 

Comme il est montre ci-dessus les statistiques C* *;^^ , pour 1 <. 
e, f, g, h, k< N sont des fonctions a six entrees, il alors possible de les 
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ranger dans > une matrice {N ^ x N ^) que Ton nommera matrice 
d'hexacovariance ffx. 

Principe mis en oeuvre dans le precede selon Tinvention 

5 Le procede selon Tinvention utilise notamment une propriete de 

multilinearite des cumulants et la gaussianite du bruit qui se traduisent 
par Texpression matricielle suivante 

Cm.x = [il®t^- 1) ® A*] Cm,s [A®^^- 1) ® il*]H (10) 

ou G/^x et Cm,s sont les matrices des statistiques d*ordre m, definies 
10 precedemment, de taille respective (iV<? x N^) et (P^ x P^), et associees aux 
vecteurs x(t) et s(f)ou A^(^ - ^) correspond a une notation adoptee definie 
airisi : la matrice B®^ designe la matrice B elevee a la puissance (au sens du 
produit de Kronecker) k, c'est-a-dire en effectuant le produit de Kronecker 
jBf^fc = B0B®,,,®B , en posant B^o = 1. 

k fois 

15 Rappelons des lors quelle est la definition du produit de Kronecker : soient* 
il et B deux matrices de taille respective (La x Ca] et (Lb x Cb) , le produit de 
kronecker D^^ A® B est une matrice de taille {La Lb x Ca Cb] definie par D = 
{Aij B)l<i< LA, l<j< CA, 

Sans sortir du cadre de I'invention d'autres ecritures associees 

20 d*autres rangemeiits des cumulants peuvent etre utilisees comiiie par 
exemple : 

. Cm.x = ] Cm,s. I [il®^ ] " ( 10 bis) 

ou 1 est choisi arbitrairement entre 1 et q et ou Cm>, i. est la matrioe des 
statistiques d'ordre m = 2q de s(t) associee a Tindice 1 choisi.* thaque 
25 ecriture conditionne le nombre de sources potentiellemerit identifiables a 
partir d'un reseau donne. 

Dans la suite de la description, I'analyse est donnee en utilisant 
Tecriture de la relation (10). 
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Dans la mesure ou les sources sont independantes, la matrice 
des statistiques d'ordre m associee aux sources, Cm.s , est une matrice 
diagonale. Neanmoins, elle s'avere ne pas etre de rang plein. Le precede 
salon I'invention considere une matrice determinee a partir d'une matrice 
diagonale de rang plein des autocumulants et de la matrice representant la 
juxtaposition des P vecteurs colonnes relatifs aux vecteurs directeurs des 
sources : 

Cm,K = A,^„,,sA^ (11) 

OU (;m.s = diag{[C/ /^' est la matrice diagonale de rang plein des 

autocumulants C;^^ d'ordre m = 2q des P sources, de taille (Px f), et ou 
A, = [ai®(.-i) ® af ... ap®W-') ® a/Sr ], de taille (iV^ x P) et supposee de rang 
plein, represente la juxtaposition des P vecteurs colonnes [ a^®(9-i) (8> ]. 
De plus, nous supposons que la matrice = [ ai® (9-2) ® ai* ... a^(g-2) ® 
ap ], de taiHe (TVfe-D x est egalement de rang plein. 

Le precede selon I'invention permet avantageusement d 'exploiter 
et d'extraire par exemple la totaUte de I'information propre aux vecteurs 
directeurs ctp des sources, contenue de maniere redondante dans la matrice 
des statistiques circulaires d'ordre m = 2q du vecteur d'observations x{t). 
Cm,x et plus particulierement dans la matrice Aq. 

Le precede comporte par exemple les etapes decrites ci-apres. L 
echantillons du vecteur x(t) sont supposes avoir ete observes et la matrice 
C„,x est supposee avoir ete estimee a partir de ces echantillons. 

Etape 1 : decomposer en valeurs singulieres la matrice Cn.x 

Calculer la racine carree C/..^^ de C,n,^ de rang plein, par exemple 
25 en decomposant en elements propres la matrice hermitienne C„,^ = Es 
Bs ou Ls et Es sont respectivement la matrice diagonale des P plus grandes 
{en valeur absolue) valeurs propres reelles de C,„,. et la matrice des vecteurs 
propres orthonormes associes. Cette etape fait apparaltre le Uen existant 
entre Cj^ et Ag, : 

30 Cnl/x^ =Es\Ls\ 1/2 = A, i;^,^ VH = [ a®„-„ ® ... ap'i'i-^) ® ap ] ^i// V» (12) 
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ou V est une matrice unitaire, de taille [P x Pj, unique pour Ls at JB, 
donnees, et ou \Ls\^^^, Cm,s /^^ sent des racines carrees respectivement de 
\Ls\ et ^m,s ( I . I designe Toperateur valeur absolue). 

Ls et Bs sont par exemple respectivement la matrice diagonale des P plus 
5 grandes (en valeur absolue) valeurs propres reelles de Cm,x et la matrice des 
vecteurs propres orthonormes associes. 

Pour une matrice A, de rang plain, il est possible de verifier que 
llij^othese (H4} equivaut a supposer les elements diagonaux de Ls non 
nuls et de meme signe. 

10 Etape 2 

Extraire de la matrice C^^i= [FiT, Tiv^JT les N blocs matrideis 
Tn : chaque bloc Tn de taille (N (9 - ^) x est constitue des iV - D lignes 
successives de CmJ a compter de la « N i^~ (rL-l)+l »-ieme. 

Etape 3 

15 Construire les N(N-1) matrices 0ni,«2 definies, pour tous 1 < m 

712 < N, par ©ni.ri2 = Tni rri2 OU # designe Toperateur de pseudoinversion. 

En notant pour tout 1 < n < iV, <E>n = diag([Ani. An2 Anp]) ou Ay est la 

composante de A situee a la i-ieme ligne et j-ieme colonne, on a Tegalite Fn 
A(qA) On CmJ pour tout 1 < n < iV, ct le fait que la matrice V diagonalise 

ft 1/2 1 1/2 

20 ' conjointement les N(N-l) matrices 0ni.n2 = Fni r„2 = V ^t^s Oni <I>„2 ^m.« V", 
qui rappelons-le, sont de taille (Px f). 

Etape 4 

Determiner la matrice Vsou solution au probleme de 
diagorialisation conjointe des N(N-l) matrices 0ni,n2 par exemple eh utilisant 
25 une methode de diagonalisation decrite dans la reference [2). La matrice 
Cm,x j^^^V^oh ou Vsoi = V T est une matrice unitaire diagonalisant 
conjointement les matrices 0rti.n2 a une matrice T triviale unitaire pres est 
une estimee de la matrice Aq a une matrice triviale pres. 
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Differentes methodes connues de 1 Homme du metier permettent 
d'extraire de Aq une estimee A de la matrice de melange A. 
Etape 5 
Etape 5A 

5 Une maniere de proceder consiste par exemple a moyenner les K 

= Ni^ - 1) blocs (S,)* de taille (N k P\ (pour tout 1 £ /c ^ iV«?- D (le bloc S^t est 
constitue des i\riignes successives de ^ = (EiT ^^t^t a compter dclsL^N 
(/c-l)+l »-ieme). ou bien de n'en retenir qu'un seul, par exemple (Z,)*. Cette 
approche permet d'estimer, dans le desordre et a une amplitude pres, les P 
10 vecteurs directeurs Op et done, a une matrice triviale pres, la matrice de 
melange A. 

Etape 5B 

Une autre maniere consiste par exemple, pour chacune des P 
colonnes bpde^= ( ai®(9-^) ® af ... a/(9-Ji 0 ap ], 

• extraire les K = vecteurs b,(k) empiles les uns en dessous des 

autres tels que : 
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. convertir les dits vecteurs colonne b,ik) = {A,,... Aj,) [a^® a,*] de taille 

X 1) en matrice B^k) = (Ap... Aj,) [a,a, h] (ou 1 < iJ<N) de taiHe [Nx 
20 N) et 

. decomposer conjointement ces K = M<^2i matrices en valeurs singulieres 
. (SVD) : le vecteur propre commun aux K matrices Bp{k) et associe a la 
plus forte (en module) valeur propre. est done un vecteur colonne de la 
matrice A. A noter que la quantite {A^... Aj,) est dans le cas present le 
25 produit de {q - 2) composantes de A. 

Cette etape de traitement sur les P colonnes bp de A^, permet 
d'estimer, dans le desordre et a une phase pres, les P vecteurs directeurs a, 
et done, a une matrice triviale unitaire pres, la matrice de melange A. 
Etape 6 
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La matrice de melange A representative des vecteurs directeurs 
des sources contient a elle seule Tinforraation necessaire a la localisation 
angulaire des sources. Dans ce contexte, a partir de Testimee des differentes 
colonnes de A, il est possible de mettre en oeuvre une methode arbitraire de 
5 goniometrie exploitant cette information. Une telle methode est par exemple 
presentee dans le document [18]. 

Etape 7 

Pour estimer le vecteur de sources s(t), en contexte surdetermine 
(i.e. lorsque P < N), le procede applique aux observations x{t) un filtre de 
10 type LIT utilisant explicitement I'estimee de la matrice de melange A, On 
peut par exemple choisir le filtre FAS decrit dans la reference [17], optimal 
en presence de sources decorrelees. 

Le procede comporte par exemple une etape 0, qui consist^ a 
15 construire, a partir des differents vecteurs d'obsej-vations x{f), une estimee 
t:m,x de la matrice de statistiques Cm^ des observations, selon la methode 
donnee precedemment. Dans ce cas les etapes 1 a 6 du procede sont mises 
en oeuvre sur Testimee tim,x de la matrice. 

20 Critere de performance 

Selon une variante de realisation, le procede comporte une etape 
utilisant un nouveau critere de performances pour revaluation 
d'identification aveugle de melange. Ce critere n'est pas global et permet/ 
d'evaluer la qualite d'identification de chaque vecteur directeur estime : on 
25 peut alors comparer deux methodes distinctes d'identification aveugle vis a 
vis de; chaque vecteur directeur, et done de chaque source. Ce qritere est 
Textension a Tidentification aveugle de melange du critere base sur le SINR 
(rapport signal sur bruit) donne dans la reference [17] introduit pour 
Textraction aveugle de sources. C'est un P-uplet decrit par 

30 D(A, A) = (ai, as, ap) (15) ou 
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«p = min [d(fl^, a.)] (16) 

et ou d(u,i^ est la pseudodistance entre les vecteurs u et v, defrnie par 

d(«,t^ = i-i(«.v)fiHi"iHr (17) 

Notons que < . , . > designe le produit scalaire defini pour deux vecteurs de 
5 meme dimension. 

Le procede decrit pour I'ordre m = 2q s'appHque notamment pour 
les statistiques a I'ordre 4 et les statistiques a I'ordre 6 par exemple selon 
les exemples donnes ci-apres. 

10 AppHcation de la methode pour la separation aveugle de sources a 
I'ordre 4 

Une variante de realisation du procede denommee ICAR 
(Independent Component Analysis using Redundancies in the quadri- 
covariance) exploite les statistiques d'ordre m = 4 (q = 2), correspondant a 
15 la matrice de quadricovariance circulaire des observations Cm,x, notee Q^. 
Ce procede permet I'identification aveugle du melange instantane A ou 
I'extraction autodidacte des sources dans le cas ou N> P, autrement dit 
uniquement lorsque le melange est surdetermine. 

Le modele (1) est suppose verifie ainsi que les hypotheses Hi^ a 

20 I'ordre 4. 

Le procede denomme ICAR exploite notamment I'expression (H) qui 
traduite pour des statistiques d'ordre 4 s'exprime : 

Qx = A2Ci.sA2» (18) 
oii C4.« = diag([Ci\i%,...,C^Pj) est la matrice de rang plein des 
25 autocumulants C^S,, d'ordre 4 des sources, de taiUe (Px f), et ou ila = [ ai i8> 
ai* ... ap® ap ], de taille (iV2 x ^) et supposee de rang plein, represente la 
juxtaposition des P vecteurs colonnes ( Up ® a$ ]. De plus, en supposant 
que la matrice de melange A, de taille (iVx Fj, est egalement de rang plein. 
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Le precede effectue par example les etapes 0 a 5 decrites dans le 
cas de Tapplication a Tordre m=2q, en utilisant les parametres suivants : 
Cm,x = et ^m.s =^4,s 

Dans cet exemple de mise en oeuvre, le precede peut aussi 
5 comprendre une etape 0 qui consiste a : construire, a partir des differents 
vecteurs dbbservations x(t), une estimee la matrice de quadricovariance 
Qx des observations. Les etapes 1 a 6 sont ensuite deroulees sur cette 
valeur estimee. 

Ebcemples de resultats obtenus en appliquant la methode a I'ordre 4 

10 Les figures 2, 3 et 4 representent dans un diagramme ou 

Tabscisse correspond au nombre d'echantLllons L et Torddnnee aux 
performances, des resultats de simulation presentant les performances du 
precede selpn Tinvention a Tordre 4 presente ci-dessus dlCAR (en mettant 
en oeuvre Tetape 5B), ICAR2 (en mettant en" oeuvre I'etape 5A) et dies 

15 methodes de separation aveugle de sources (COMl, COM2, JADE, FastICA) 
connues de lHonmie du metier. Les conditions de simulation sont les 
suivantes : 

• P = 3 sources non filtrees, une BPSK et deux QPSK, sont supposees etare 

e. 

receptionnees sur un reseau circulaire de N = S capteurs tel que R/X = 
•20 0,55 (avec R et X le rayon du reseau et la longueur d'onde) et tel que le 
rapport signal sur bruit, note SNR, soit egal a 20 dB pour chaque 
source. ' ■ . ' 

• Le bruit est Gaussien et non correle spatialement. 

^ • Les trois sources sont en bande de base et leur temps symbole est choisi 
25 ^ egal au temps d'echantillonnage. 

• Le critere employe pour apprecier au mieux les resultats d'extraction de 
la source p pour une methode donnee est le rapport signal sur bruit plus 
brouilleur maximal associe a la source p, plus connu sous le 
denominatif de SINRMp [17].I1 peut alors etre compare au SINRMp 

30 optimal calcule en utilisant non pas la matrice de melange estimee mais 
au contraire la matrice de melange exacte ainsi que les statistiques. 
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exactes des observations. C'est cette comparaison qui est presentee dans 
les figures 2-4. 

Plus particulierement, la figure 2 represente les SINRM de la 
source 1 associes a ICAR, ICAR2 mais aussi aux methodes de separation 
aveugle de sources les plus performantes actuellement telles que JADE, 
COMl, COM2 et FastlCA. 

La figure 3 represente le SINRM de la source 2 pour les memes 
methodes (ICAR, ICAR2, JADE, COMl, COM2, FAST ICA)et la figure 4 ceux 
de la source 3. 

Sur chacune des trois figures, apparait le fait que les deux 
methodes ICAR et ICAR2 ont des performances tres bonnes et legerement 
meiUeures que les autres methodes JADE, COMl, COM2 et FastlCA. Quant 
a I'algorithme FastlCA, sa meilleure performance conceme la source 2 pour 
laquelle il converge completement a partir de 550 echantillons. 

Quant a la difference de resultat entre ICAR et ICAR2 elle s'avere 
dans cette configuration etre negUgeable devant I'ecart de performance entre 
les methodes ICAR et les methodes JADE, COMl, COM2 et FastlCA, 
toutefois tres bonnes. 



Application du precede k Tordre 6 

Selon une autre variante de realisation, le precede utilise les 
statistiques a I'ordre 6. Cette variante denommee BIRTH (Blind Identification 
of mixtures of sources using Redundancies in the dcCTa Hexacovariance 
matripdi permet a partir d'un reseau de N capteurs d'identifier les vecteurs 
directeurs d'au plus P= sources (cas d'un reseau a capteurs differents). 
Elle permet egalement I'extraction d'au plus P= N sources dont les vecteurs 
directeurs sont explicitement identifies. 

Le precede BIRTH utilise certaines statistiques d'ordre 5, stockees dans la 
matrice dTiexacovariance C6,x designee . Ainsi, cette variante de mise en 
oeuvre exploite pleinement Tinformation propre au melange instantane A 
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des sources, Gontenue dans Hx, notamment grace a une astucieuse ecriture 
de Hx par rapport aux vecteurs directeurs des sources, et ce au moyen de la 
propriete de multilinearite des cumulants : 

Hx = i43 C6,sA3H (20) 

ou i^6.s = diag([Ci;^^^^^^^^ est la matrice de rang plein des 

autocumulants C^$$,s d'ordre 6 des sources, de taille (Px P), et ou ila = [ 
ai ^ ,®2® ai ,* ap ,®2® ], de taHle (JV^ x et supposee de rang 
plein, represente la juxtaposition des P vecteurs colormes Ic^® ^ <g) Op * ] = 
[dp® €Zp ® aj? ]. De plus, nous supposons que la matrice A2 = [ ai ® ai* 
ap <S) ap ], de taille (iV^ x f), est egalement de rang plein. 

Le precede a Tordre 6 comporte les etapes 1 a 6 et Tetape 0 
decrites ci-dessus en utilisant les parametres suivants : Cm,x ^ Hx et in-3. 

Simulations 

Les figures 5, 6, 7 et 8 representent dans un diagramme ou 
Tabscisse correspond au nombre d*echantillons L et Tordonnee aux 
performances, le bon fonctionnement du precede selon Tinvention a Tordre 
6. Les conditions d'obtention des courbes sont les suivantes : 

• supposons que P = 3 sources statistiquement independantes, plus 
particulierement deux QPSK et une BPSK toutes trois nbn filtrees, sont ' 
regues sur un reseau lineaire de iV = 2 capteurs tel que R f X ^ 0.55 
(avec R et X etant respectivement le rayon du reseau et la longueur 
d'onde). , 

Les trois sources, supposees synchronisees, ont meme rapport signal 
sur bruit/ note SNR et egal a 20 dB pour chaque source, avec pour tei^ps 
s)rmbole quatre fois le temps d'echantillonnage. v 

• La BPSK est choisie en bande de base, alors que les.deux QPSK ont des 
porteuses respectivement egales a la moitie et au tiers de la frequence 
d'ech antillonnage . 
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• La matrice de melange A est choisie de maniere a ce que les vecteurs 
colonnes de la matrice soient lineairement independants. Le bruit est 
quant a lui Gaussien et non correle spatialement. 

• Le melange instantane de sources bruitees est considere comme un 
5 melange surdetermine car le nombre de sources est superieur au 

nombre de capteurs. Les algorithmes JADE, COMl, COM2, S3C2 
connus de momme du metier et le procede selon invention a I'ordre 6 
sont mis en oeuvre afin d'identifier de maniere autodidacte le melange 
surdetermine A. 

10 • Le critere de performances Op, defini par I'equation (18), est calcule sur 
200 realisations et ce pour chaque source p (1 ^ p < 3) : il va ainsi 
permettre la comparaison des cinq methodes. 

Sous les hypotheses precedentes, la figure 5 montre les 
variations de la quantite as a la sortie des algorithmes JADE, COMl, COM2, 
15 S3C2 et le procede selon I'invention BIRTH en fonction du nombre 
d'echantiUons. Le procede selon I'invention permet d'identifier le vecteur 
directeur en question contrairement aux methodes de I'art anterieur. 

La figure 6 montre, dans le meme contexte, les variations du 
triplet D[A, A) = (ai, a2, as), associe au procede selon I'invention a I'ordre 6, 
20 en fonction du nombre d'echantiHons. Les trois coefficients ap decroissent 
rapidement vers zero a mesure que le nombre d'echantillons augmente. 

La figure 7 presente les variations de la quantite as a la sortie des 
methodes selon I'art anterieur JADE, COMl, COM2, S3C2 et le procede 
selon I'invention en fonction, cette fois, du rapport signal sur bruit (SNR) de 
25 la source 3. La methode BIRTH, parvient parfaitement a identifier le vecteur 
directeur de la source 3 meme pour une faible valeur de SNR. 

Finalement, supposons que les P = 3 sources precedentes soient re?ues sur 
un reseau circulaire de iV= 3 capteurs tel que RfX = 0.55. 
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La figure 8 montre alors les variations de la quantite az a la sortie 
des algorithmes JADE, COMl, COM2 et BIRTH en fonction du nombre 
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d'echantillons- : Talgorithme BIRTH fonctionne egalement en contexte 
surdetermine, c'est-a-dire lorsque le nombre de sources est inferieur au 
nombre de capteurs, et bien que des cumulants d'ordre 6 des observations 
doivent etre estimes, la vitesse de convergence de BIRTH est du meme ordre 
de grandeur que celle des methodes citees plus haut. 
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1 - Precede d Identification aveugle de sources au sein d'un systeme 
comportant P sources et N capteurs caracterise en ce qu'il comporte au 
moins une etape d Identification de la matrice des vecteurs directeurs des 
sources a partir de lUnformation propre aux vecteurs directeurs ap des 
sources contenue de maniere redondante dans les statistiques circulaires 
d'ordre in=2q du vecteur des observations re?ues par les N capteurs. 



2 - Procede selon la revendication 1 caracterise en ce que les statistiques 
circulaires dbrdre m - 2q s'expriment en fonction d'une matrice diagonale 
de rang plein des autocumulants des sources et d'une matrice representant 
la juxtaposition des vecteurs directeurs des sources de la maniere 
15 suivante ; 

Cm,^^ AqCm,sAq^ (11) 

ou Cm,s = d^^g{[Cl;l -j^^ >--'Cp'pZ'.p^^ ]) ^st la matrice diagonale de rang plein des 
autocumulants C^^:^ dbrdre m - 2q des sources, de taille (Px Pj, et oii 
= [ ai®(9-i) (8) ai* ... a/(^i) ® a? ], de taiUe (iV^ x et supposee de rang 
20 plein, represente la juxtaposition des P vecteurs colonnes [ a^®(g-i) ® ], 

3 - Procede selon Tune des revendications 1 et 2 caracterise en ce qu'il 
comporte au moins les etapes suivantes : 

Q : construire, a partir des differents vecteurs d 'observations x(tj, une 
25 estimee C^^ de la matrice de statistiques Cm^ des observations, 

i : decomposer en valeurs singuUeres la matrice C^^ , en deduire une 
estimee P du nombre de sources P et une racine carree C^'^ de C n«r 



I CI UCjJUlJ 
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exemple en prenant C^'^, = JSs \Ls\^^^ ou |.| designe Foperateur valeur 
absolue, ou Ls et Es sont respectivement la matrice diagonale des P plus 
grandes (en valeur absolue) valeurs propres reelles de C^^ et la matrice des 
vecteurs propres orthondmies^^^a^^^ 

5 2 : extrake^ la matricexNCJJJ^^ =1 ^(T i"^ les N blSiC^matriciels Tn : 



cha^ue bloc<r>it' dd^^^ (iV (9 - i) x f) est constitue des JV lignes 
sticcessiyesvde a compter^drla «^ - il(rt- 1 ) + 1 »-ieme iv/ 




// 3 :/eonstnaire les ^Mi^^^ matrices ,©ni.r.2;definies,vpoui%|ous 1 <-n4^^m <^^\ 

if . /// ;i r- \\\ \ 

/ 10 / 4 Txdeterminer la matrice ysoi, solution' aus probleme dfeXdiagonalisatibn 

; / / / wy'^'-'i-/. ^' \^ -. \\ \ ( I 

Vconjointe des iV^[^^ \\ \ — 



5: pour /chacune des P colorliies. bp /de^^^^^^ K= .N(^^) vecteurSt 

yZ""'- ^ ' ' ^ \ ^ L ^ ~ 7 ^ 1 H ^^^^ 

bp(/c) empiles les uns en dessous des^autres dans le vecteur jftpf = [bp{l)'^,2i 



15 A 6 : cdnvertir les^^^ vecteurs colonne i&pffc) de taille [N^ x 1) en niatrice 
2\: DecomposerU^ en^aleurs sirigulier:es,^o^^ 




\ X.^r'- 



conjointement les\ff iW^^r matrices Bp(/c) en et retenir comma estimee 



des^ecteurs colonnes de A, le vecteur propre commun^auxK matinees 



20\ asspjcie^a la plus forte (en module) valeur propre 



8\ repeter les etapes 5 a T^pour— chacune des P colonnes de 
d'estimer, dans^le^esordre et a une phase pres, les^P vecteurs directeurs cu 
et donc>aHine matnce tnviale^unitair^ pres^ \la matrice de melange A.. 




25 4 ~ Procede selon Tune des revendications 1 a 3 caracterise en ce que le 
nombre de capteurs N est superieur ou egal au nombre de sources P et en 
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ce qu'il comporte une etape d'extraction des sources consistant a appliquer 
aux observations x{t) un filtre construit a I'aide de I'estimee AdeA. 



5 - Procede selon IHine des revendications 2 a 4 caracterise en ce que C,^ 
5 est egale a la matrice de quadricovariance Qx et en ce que m = 4. 



6 - Procede selon l\ine des revendications 2 a 4 caracterise en ce que 
est egale a la matrice dliexacovariance Hx et en ce que m = 6. 

0 7- Procede selon IHine des revendications 1 a 6 caracterise en ce qu'il 
comporte une etape d'evaluation de la qualite de I'identification du vecteur 
directeur associe en utilisant un critere tel que 

D(A,A) = (ai, as, ap) 

ou 

et ou d{u,tij est la pseudo-distance entre les vecteurs u et v, tel que : 
d{u, = 1 - 



u v\ 



8 - UtiUsation du procede selon I'une des revendications 1 a 7 a un reseau 
de communication. 



9 - UtiUsation du procede selon I'une des revendications 1 a 7 pour de la 
goniometrie a partir des vecteurs directeurs identifies. 
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